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1. INTRODUCTION 


This report is the Final Technical Report on NASA Grant NAGW-240, 
monitored by Drs. Clint Brown, Gary Hicks and Randolph Graves. The report 
describes experimental work on turbulent, supersonic shear layers performed 
at the Gas Dynamics Laboratory of Princeton University during the period 
September 1, 1981 to August 31, 1984. 

The behavior of turbulent shear flows at high Reynolds numbers and 
supersonic speeds is of great practical interest, particularly in the internal 
and external aerodynamic design of aircraft. A good example is given by the 
flow through the compressor and turbine. Here, very complicated flows occur, 
and the boundary layers can experience severe adverse pressure gradients, 
interactions with shock waves, longitudinal curvature and possibly separation 
and reattachment. Similar observations can be made regarding the disturbances 
suffered by the external flow over the fuselage, wing and control surfaces. 

Unfortunately, our present understanding of these complex flows is 
rather limited. As a result, the computations display severe shortcomings, 
and these inadequacies were well documented at the recent Stanford Conference 
on Complex Turbulent Flows (Kline, et.al. 1981). Consider the example shown 
in Fig. 1. Here, the calculations by Visbal and Knight (1983) are compared 
with the measurements by Settles, et.al. (1977) in four two-dimensional 
compression corner flows at Mach 3. The flows range from attached (8°), to 
incipient separation (16°) to separated (20°, 24°). The computer code was 
based on the Reynolds-averaged Navier-Stokes equations using the Baldwin-Lomax 
algebraic eddy viscosity model. As can be seen, the comparisons are not very 
favorable, especially for the higher corner angles and the results clearly 
demonstrate the inadequacy of the turbulence model. 
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For three-dimensional interactions, similar calculations have resulted 
in rather better agreement with experiment, at least in some cases. For 
example, Knight (1983) found good agreement between his calculations and the 
experimental study of a sharp fin interaction performed by Oskam (1976). Such 
calculations, however, require careful validation before they can be considered 
reliable and accurate. The recent computations by Horstman (1984) demonstrate 
how cautiously we must proceed. Horstman calculated a large number of swept 
compression corners, systematically varying the corner angle a and the sweep- 
back angle A. These results were compared with the experimental work by Teng 
and Settles (1982), Settles and Bogdonoff (1982) and Settles, et al. (1984), 
and Horstman found good agreement over a considerable range of a and A. At 
high sweepback angles, however, sharp discrepancies occurred. In addition, 
the computations missed some important flow details at lower angles. Thus, 
computations can appear to be successful over some range of parameters and 
yet display severe shortcomings outside this range; without appropriate 
experimental work such discrepancies might pass unnoticed. 

Experimental work, therefore, is required (i) to understand the physical 
mechanisms which determine the behavior of the flow, (ii) to develop improved 
turbulence models, and (iii) to generate data sets for the validation of compu- 
tational work. In particular, we require further turbulence data; the existing 
measurements are scarce and generally of dubious quality (see, for example, 
Fernholz and Finley, 1981). 

These needs were recognized by NASA, and the Gas Dynamics Lab was 
asked to begin a long range experimental program to study the behavior of 
turbulence in supersonic shear flows. We feel that this program has achieved 
major progress in the understanding and documentation of supersonic shear 



layer behavior, and the current report summarizes our achievements. Further 
details are available from the publications produced under NASA sponsorship, 
and these are listed in Appendix A. 

The research can be conveniently divided into three major areas: 
development and improvement of turbulence measuring ;echniques for supersonic 
flows (see Section 2); detailed experimental investigations of a limited number 
of flow geometries (Section 3), and development of physical models to explain 
the observed behavior of the flows (Section 4). In addition, a number of other 
activities took place which are not so easily classifiable, and these activities 
are described in Section 5. The conclusions are given in Section 6, together 


with some recommendations for future work. 



2. INSTRUMENTATION DEVELOPMENT AND TESTING 


Previous turbulence data displayed considerable scatter (see, for 
example, Fernholz & Finley, 1981), and a major effort was therefore under- 
taken to improve the accuracy of these measurements. Two techniques are 
available: hot-wire anemometry and laser-Doppler velocimetry. Early studies 

using the LDV technique (which requires seeding of the flow) showed that this 
method suffered from lack of resolution and seeding non-uniformity close to 
the wall and in strong velocity gradients. There was considerable doubt 
whether the particles followed the very sharp changes in flow direction which 
occur through such interactions. Our previous experience had shown that the 
hot-wire, with some restrictions, could be used for examining these flows, and 
our attention was therefore focussed on improving hot-wire techniques, and 
developing them to the point where hot-wires could be used to their full capa- 
bility under the harsh conditions of supersonic tests. 

In a blowdown facility such as the 8" x 8" supersonic channel at the 
Gas Dynamics Laboratory, run times are typically limited to one or two minutes. 
Such short run times make the use of conventional constant current hot-wire 
anemometers impractical. These anemometers require careful adjustment of the 
frequency response under actual operating conditions, and data acquisition is 
therefore rather slow. Thus, we have concentrated our efforts on improving 
constant temprature anemometer techniques. These systems compensate for fre- 
quency response automatically and considerably improve data acquisition effi- 
ciency. We have succeeded in developing reliable calibration and data 
acquisition techniques for normal and inclined wires operated in the constant 
temperature mode, and these methods have been reported at several conferences 
and in a number of archival journals (see Appendix A). To summarize, significant 



original contributions in this area include the design of reliable and 
durable hot-wire probes, the analysis and experimental validation of inclined 
wire sensitivity in supersonic flow and the identification of the effects of 
normal Mach number, end conduction and wire bowing on inclined wire sensitivity. 

| It was also shown, for the first time, that constant temperature systems dis- 
play significant non-linearities at low overheat ratios, and they are therefore 
unsuitable for the measurement of temperature fluctuations. Only mass-flow 
fluctuatio- s, (ou)', can be satisfactorily measured with a constant-temperature 
system . 

Two oifferent probes were designed and tested: a normal wire probe, 

for measuring (pu) ' , and an inclined wire probe, for measuring the mass-weighted 
Reynolds shear stress (pu)'v'. The probes, and the calibration and testing 
procedures, are described briefly below. 

2.1 Normal Wire Probe Design, Calibration and Testing 

The current probe design is the result of considerable development. 

In our first design, bare tungsten wire was spot-welded to the prongs using 
a tungsten electrode. This proved to be highly unsatisfactory. Not only was 
the active wire length subjected to aerodynamic interference from the relatively 
bulky prongs, but it was also difficult to achieve a satisfactory bond between 
wire and prong, and wire breakages were very frequemt. Instead, our current 
probe design closely follows that recommended by Kovasznay (1950). The tungsten 
wire (5 ym diameter) is first electroplated with copper and then soft-soldered 
to the prongs. A dilute suphuric acid solution is then used to etch away the 
copper coating and expose an active portion of tungsten wire approximately 
0.8 mm long (Fig. 2). To avoid strain-gaging, a small amount of slack is 
usually introduced. This probe design drastically reduces wire breakage, and 
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in addition, minimizes the interference effect of the bow shocks emanating 
from the tips of the prongs (Figs. 3 and 4). 

The probe was connected to a DISA 55M10 constant temperature anemometer. 
The overheat ratio was varied by changing the bridge resistance. The frequency 
response (deduced from, a square-wave test) was optimized by adjusting the ane- 
mometer gain and filter setting. All wires were checked for strain-gaging, 
and those found to be suspect were discarded. 

The anemometer output was separated into a mean and a fluctuating 
component by low- and high-pass filters each set at 10 Hz. The fluctuating 
component was digitized directly at 500 kHz sampling rate by a Preston Scientific 
GMAD-1 A/D converter, and the raw data was stored on-line in the memory of a 
Hewlett-Packard HP1000 minicomputer for further processing. The mean component 
of the output voltage was also recorded, along with other mean quantities, by 
a slower A/D converter. 

The wires were tested and calibrated in a small Mach 3 pilot tunnel 

with a working section measuring 49.3 mm x 74.5 mm. The stagnation pressure 

5 5 2 

was varied between 4 x 10 and 14 x 10 N/m which, for a 5 pm wire, gave a 
wire Reynolds number range of approximately 80 to 250. 

We demonstrated that the constant-temperature hot-wire anemometer is 
inherently unsuitable for measuring turbulent temperature correlations; the 
major reason is the non-linearity of the temperature sensitivity at low 
overheat ratios. The instrument is therefore restricted to measurements of 
the mass-flow fluctuations. If temperature fluctuations are present, high 
overheat ratios are desirable to avoid contamination of the mass-flow signal 
by contributions from the fluctuating temperature. Very high overheat ratios 
may be required if we wish to ignore these contributions entirely. 
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The maximum frequency response of the system depends on the anemometer 
roll-off frequency, the spatial resolution of the probe and the maximum A/D 
conversion rate. The maximum frequency response required depends on the 
experiment; all requirements on the frequency response become less stringent 
as the typical size of the shear layer increases. 

It was found that the static calibration of the anemometer could be 
adequately represented by a modified King's Law. This calibration is a 
function of mean stagnation temperature and corrections are required to 
account for this dependence if the stagnation temperature varies with time, 
or with position in the flow field. The corrections were found to be signifi- 
cant, and a satisfactory correction procedure was suggested. 

To demonstrate the constant-temperature hot-wire technique in practice, 
some measurements were made in the boundary layer developing on the tunnel 
floor of the Princeton University 203 mm x 203 mm Supersonic Wind Tunnel. 

The freestream Mach number was 2.9, the wall conditions were near-adiabatic, 

and the rms mass-flow turbulence level in the freestream was approximately 1%. 

5 2 

The tunnel was operated at a stagnation pressure of 6.9 x 10 N/m , which 
gave a unit Reynolds number of 6.3 x 10^/m. At the measuring position, the 
boundary layer thickness was about 26 mm, with a Reynolds number based on 
momentum thickness of 77,600. 

Consider the data presented in Fig. 5. Figure 5a shows the mass-flow 
turbulence intensity as measured by three different wires, each operated at 
an overheat ratio of approximately 1.0. The repeatability of the measurements 
is obviously very satisfactory. The effect of varying the overheat ratio is 
demonstrated in Fig. 5b. Clearly, the data reach an asymptotic level as the 
overheat ratio increases, suggesting that the results taken at high overheat 
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are not significantly "contaminated" by total temperature fluctuations. 

The inferred velocity fluctuation intensity is shown in Fig. 5c. 

The present measurements appear to be a little higher than most comparable 
data, particularly near the wall. The relatively high levels shown by our 
results may be due to the good spatial and temporal resolution of our measure- 
ments; significantly, our measurements agree well with the data of Johnson and 
Rose (1975) who used both LDV and hot-wire systems. 

In conclusion, the overall trend shown by our data, and the level of 
quantitative agreement with previous data appears to be satisfactory. 

2 .2 Inclined Wire Probe Design, Calibration and Testing 

The inclined wire probe design was similar to that used for the normal 
wire probe described in Section 2.1. The tungsten wire (5 pm diameter) was 
first electroplated with copper and then soft-soldered to the prongs. The 
central portion of the copper coating was etched away to expose an active 
portion of the tungsten wire approximately 0.8 mm long (see Fig. 6). For an 
inclined wire it is particularly important to isolate the central active 
length from the flowfield interference caused by the prongs, and the copper 
plated stubs served this purpose admirably (see Fig. 7) . The angle formed 
between the mean flow direction and the normal to the wire for our probes 
was either 30°, or 45° (approximately). 

The calibration facility, anemometer equipment and data analysis 
techniques were as described in Section 2.1. The inclined wire was calibrated 
for mass-flow rate sensitivity according to the procedure used for a normal 
wire. To calibrate the inclined wire for directional sensitivity, the probe 
was yawed through an angle of about ± 10° from its null position using the 
device shown in Fig. 8. This device yaws the wire without changing its position 


in the flow. 
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We established a reliable calibration procedure for inclined wires 

operating in supersonic flow. The longitudinal mass-flow sensitivity at a 

fixed yaw angle was found by fitting a modified King's Law to the data. 

This data correlation was identical to that used in the analysis of normal 

* 

wire data, with both correlations using an exponent n = 0.55. Several 
observations were made regarding the relative sensitivity £, defined as the 
ratio of the transverse to longitudinal mass-flow sensitivity (see Fig. 9). 
First, £ appears to be independent of the Reynolds number. Second, although 
£ depends strongly on the flow direction relative to the wire, the functional 
dependence is not simple; for instance, the cosine cooling law does not hold 
in supersonic flow. It seems best to find £ by direct calibration. Third, 
when the relative sensitivity is a strong function of the yaw angle, the wire 
should be discarded to avoid errors due to nonlinearities. Fourth, in a super- 
sonic flow £ takes higher values than commonly encountered in subsonic flows. 

The Mach number dependence of the inclined wire response was not inves- 
tigated. Present indications are that this Mach number dependence is small, 
in agreement with the findings of Reshotko and Beckwith (1958). Further work 
on this question is currently in progress. 

The calibration is a function of probe alignment and mean stagnation 
temperature. The corrections for temperature drift are similar to those used 
in the normal wire analysis. When the probe is misaligned by only a small 
amount, the senstivities may be corrected appropriately, but the measurements 
can not be interpreted without knowing the intensity of the longitudinal mass- 
flow rate fluctuations. 

Two major restrictions limit the usefulness of the inclined wire in 
supersonic flow. The first restriction requires that the normal Mach number 


I 


- 10 - 


t 


) 


must exceed 1.2. This places a limit on the maximum turbulence intensity 
that can be measured accurately. The second restriction is the necessity 
for a high system roll-off frequency. This lower limit on the frequency 
response is considerably higher than that required for normal wire measurements, 
and it may not be attainable in practice. 

To demonstrate the inclined wire technique in practice, the mass-weighted 
Reynolds shear stress (pu) 'v' was measured in a Mach 2.9 turbulent boundary 
layer. By comparing the results with corresponding data for similar flows, 
some indication of the measurement accuracy was obtained. 

In these measurements, two different inclined wires were used, and the 
angle calibrations of these particular wires were those shown in Fig. 9. Both 
wires had a nomi.ial diameter of 5 pm, a length-to-diameter ratio of approxi- 
mately 180, and in each case the frequency response in the freestreara was about 
125 kHz. The boundary layer flow was typical of a zero pressure gradient layer, 
and the flow was identical to that investigated using normal wires. Further 
details of the flowfield, and the data acquisition system were given in 
Section 2.1. 

The results are shown in Fig. 10. The Reynolds shear stress u'y' was 
deduced by assuming that (i) the pressure fluctuations were small, (ii) the 
total temperature gradient was small, and (iii) the turbulent Pradtl number 
was unity (Mikulla, in: Fernholz and Finley, 1981). 

The normal Mach number restriction affects the measurements near the 
wall, and it will mean that (cu)'v’ and u'v' are underestimated close to the 
wall. For a straight 30° wire in our Mach 2.9 boundary layer, the results 
are underestimated for y/6 < 0.1 and for a straight 45° wire the results are 
underestimated for y/o < 0.4. Khen the wire is bowed, the situation is. 
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considerably worse. For instance, with a bow of ± 10°, the results for a 
nominally 30° wire are underestimated in the region y/6 < 0.2, and a nominally 
45° wire cannot be used for y/6 < 0.6. Thus the results for the twp (slightly 
bowed) wires shown in Fig. 10 can be explained in terms of the normal Mach 
number criterion. 

The u'v' measurements shown in Fig. 10b may be compared with the "best 
estimate" of the supersonic shear stress distribution suggested by Sandborn 
(1974) . The agreement is rather poor for the 45° wire, but the 30° wire 
gives excellent agreement for y/6 > 0.2. In fact, considering the size of 
the error bars, the agreement is rather better than expected. The discrepancies 
displayed by the 45° wire for y/6 < 0.6, and by the 30° wires for y/6 < 0.2 may 
be explained in term’s of the normal Mach number criterion, as has already been 
discussed. 

In practice, hot wires are often bowed, either because of thermal 
expansion; or because of a deliberate attempt to reduce strain-gauging effects. 
In addition to bowing, end-conduction effects may be responsible for the scatter 
observed in the directional sensitivity (see, for example, Fig. 9c). 

The directional sensitivity of a bowed, inclined hot-wire was there- 
fore investigated using a simple model for the convective heat transfer. 

The static response was analyzed for subsonic and supersonic flows. It was 
shown that the effects of both end conduction and wire bowing are greater in 
supersonic flow. Regardless of the Mach number, however, these two phenomena 
have distinctly different effects; end conduction apepars to be responsible 
for reducing the non-linearity of the response, whereas bowing increases the 
directional sensitivity (see Fig. 11). 




It was also possible to calculate the temperature distribution along 
the bowed wire, and sample calculations are shown in Fig. 12. The effect of 
yaw on a straight wire is surprisingly small, and the temperature distribution 
is virtually unaffected by yaw angle. In contrast, the bowed wire temperature 
distribution is a strong function of yaw angle. 

The position of the temperature peak occurs along the portion of the 
wire with the greatest sweepback, as may be expected. In addition, the 
maximum temperature for a bowed wire is considerably higher than the maximum 
for the straight wire. This observation has important consequences for hot 
wire filaments made of tungsten; to avoid oxidation the temperature at any 
part along the wire should remain below 600°K. 

2.3 Current Development Work 

Current development work in the area of hot-wire anemometry includes 
a micro-schlieren flow visualization investigation of the flow in the vicinity 
of the wire (preliminary results are shown in Figs. 4 and 7), the design of 
a suitable crossed-wire probe, the study of the hot-wire response in transonic 
flows and the development of a probe with an array of wires for the purpose of 
gathering data simultaneously at a number of points in the flow. 

In addition to hot-wires, wall pressure transducers have been used in 
our lab for some time, particularly in the investigation of unsteady flows. 
Recently, we have extended the use of these transducers to multiple channels 
to measure space-time correlations, and to operate in combination with hot- 
wire probes to measure wall-pressure/velocity correlations. The validity 
of these measurements is restricted by the upper frequency response of the 
pressure transducers and the phase shift performance. These aspects are 
currently under investigation. 



3 . EXPERIMENTAL STUDIES 


The hot-wire development resulted in a clear understanding of normal 
and inclined wire behavior, and it enabled us to take high quality turbulence 
measurements in supersonic flows. Great progress was achieved in using these 
hot-wire techniques to investigate a variety of interesting flows. These 
flows may be conveniently distinguished into two groups: shock-wave/boundary 

layer interactions produced by two-dimensional compression corners, and 
"isentropic compressions" where the boundary layer flows over a concavely 
curved surface. These flows were chosen to investigate the behavior of 
supersonic turbulent boundary layers in moderate to severe pressure gradients. 

In addition, the reattachment of a boundary layer separating off a backward- 
facing step was investigated in considerable detail to improve our understanding 
of the reattachment process. We will now consider each group of flows in turn, 
beginning with the reattachment study. 

3.1 Reattaching Shear Layer Study 

Using a backward facing step geometry, the incoming boundary layer was 
allowed to develop over a flat plate at Mach 2.9 before separating off the 
step and reattaching on a 20° ramp. The ramp was designed to eliminate the 
lip shock at the point of separation, and the separated flow formed a nominally 
self-preserving free shear layer without extraneous shocks and disturbances, 
and therefore it provided a well-defined initial condition for reattachment. 

The test model is shown in Fig. 13. A turbulent boundary layer 
developed initially on the upstream flat plate (229 mm long). At the point 
where the boundary layer separated over a backward facing step, the boundary 
layer thickness 6 was about 3.5 mm with a Reynolds number based on momentum 
thickness of approximately 14,000. The free shear layer, formed by the 
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separating boundary layer bridged a 25.4 mm deep cavity before reattaching 
on a plane ramp (see Fig. 14) . 

The flowfield steadiness was investigated by Settles et.al. (1983) 
using microsecond shadowgrams. In contrast to subsonic reattachment, where 
the reattachment line can move significantly, only a slight "tremble" of 
the wave system was observed. The magnitude of this wave motion was small 
compared to the average shear layer thickness and no large scale unsteadiness 
seemed to be present. 

The upstream turbulent boundary layer separates without deflection at 
the corner of the backward facing step, forming a free shear layer. The 
shear layer mean velocity profiles achieve self-similarity at station SH43, 
which is approximately 17 6 0 downstream of the step. The growth of the shear 
layer is faster on the low-speed side than on the high-speed side, and, as a 
result, reattachment occurs at a point slightly below the geometric extension 
of the flat plate on the ramp surface. The static pressure rises before re- 
attachment and continues to rise well downstream. The compression wave asso- 
ciated with the shear layer curvature coalesce to form a shock wave in the 
freestream, at some distance outside the shear layers. 

Downstream of the mean reattachment point (station R27) a new boundary 
layer begins to develop. The adverse pressure gradient in this region gradually 
reduces and becomes negligible at about station R44. The mean velocity near the 
wall rapidly increases with downstream distance and the wall-layer thickness 
grows quickly. The velocity profiles initially display a very large wake 
component (as may be expected) but this soon decreases. By station R42, the 
mean velocity appears to dip below the standard logarithmic law, suggesting 
that the length scale near the wall is abnormally large. The relaxation of 



the boundary layer is certainly not monotonic. For instance, both the Clauser 
parameter G and the wake parameter it reach high values near reattachment but 
undershoot their respective equilibrium values of 6.8 and 0.55 downstream. 

This behavior agrees with that observed in many subsonic reattachment studies 
(Bradshaw and Wong, 1972; Chandrsuda and Bradshaw, 1981). 

The turbulence behavior is in some ways even more spectacular. Figure 15 
shows the downstream evolution of the maximum rms mass-flow fluctuation inten- 
sity <(pu) •>. In this figure , the turbulence intensities are normalized by 
the mean freestream mass flow upstream of the reattachment to show the ampli- 
fication of the absolute intensity. To interpret the results in terms of the 
local freestream mass-flow rate, it may be noted that the freestream mass-flow 
rate increases by 97% through the compression. The arrows in Figs. 3 and 4 
indicate the sonic point locations. The results below these points should 
be ignored, since the hot-wire calibration is not valid in the transonic and 
subsonic regions. 

As can be seen, the fluctuation intensity in the free shear layer 
increases slowly with downstream distance. In contrast to the mean flow 
measurements, it appears that the turbulence intensity profiles do not achieve 
similarity at any stage before reattachment. However, this observation may be 
incorrect. The relative temporal and spatial resolution of the hot wire improves 
improves with increasing shear-layer thickness, and this may give the erroneous 
impression that the turbulence intensity is increasing. 

The maximum turbulence intensity rises rapidly as the shear layer 
approaches the ramp, and at reattachment it reaches a level of almost 40% 

(see Fig. 15). Downstream of this point, the intensity continues to rise 
before reaching a maximum at station R42 (which coincides closely with the 
point where the static pressure gradient becomes negligible). 
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This large increase in turbulence level, persisting downstream of 
reattachment, is rather unexpected. In subsonic flow, downstream of a 
backward facing step, it has been observed that the turbulence level reaches 
a peak before reattachment and then decays rapidly. In contrast, in super- 
sonic flow, we observe a dramatic amplif ication of the turbulence level. 

This experiment, therefore, clearly demonstrates the effect of compressi- 
bility on the turbulence behavior. 

By using an approximate form of the turbulent kinetic energy equation, 
we demonstrated that mean dilatation significantly contributed to the 
amplification of the turbulence intensities. The overall flow behavior, 
however, was far from simple, and many competing influences were present. 

In particular, the turbulence length scale appeared to be both reduced by 
eddy bifurcation near reattachment and amplified by the action of extra 
strain rates due to dilatation and longitudinal curvature. In our analysis 
we ignored the possibility of turbulence amplification caused by unsteady 
oscillation of the wave system, and by local deformation of the compression 
waves by the turbulence itself. Several authors (Anyiwo and Bushnell, 1982; 
Zang et al., 1982) have suggested that these mechanisms may be important. 

In short, successful prediction of the present flow will require some very 
sophisticated modeling, and the challenge to the predictor is clear. 

3.2 Compression Corner Studies 

Three compression corners were surveyed using both normal and inclined 
wires. The corner angles were 8°, 16° and 20°, which, at a Mach number of 2.9, 
correspond to a fully attached flow, a flow on the point of separation and a 
separated flow with a small separation bubble. The incoming flow was identical 
in all three cases, and the upstream boundary layer <$ 0 was approximately 



26 mm thick, with a Reynolds number based on momentum thickness of 77,600. 

The results were reported in publications [l], [3], [l2], [13] and [17] 

(see Appendix A), and the data were tabulated in publications [9], [l4], [15] 
and [20]. These data compilations are particularly useful in disseminating 
results to computors (see Section 4), and we have adopted their production 
as standard practice. 

The tunnel configuration is shown in Fig. 16. The tunnel was 

5 2 

operated at a stagnation pressure of 6.9 x 10 N/m , and the unit Reynolds 
number was 6.3 x 10^/m. The wall conditions were near-adiabatic and the 
freestream rms turbulence intensity was approximately 1%. The models were 
mounted on the tunnel floor and fitted with aerodynamic fences to avoid 
interference from the side-wall boundary layers. Blockage problems limited 
the model lengths downstream of the corner to 195 mm for the 8° model, 151 mm 
for the 16° model and 121 mm for the 20° model. 

Constant temperature hot-wire anemometers were used throughout, and 
details of the operating procedure were given in Sections 2.1 and 2.2. 

The normal wire measures fluctuations in mass-flow rate (pu)'. To 
determine the behavior of the fluctuations in velocity, we invoked Morkovin's 
(1962) "Strong Reynolds Analogy," which assumes that the pressure fluctua- 
tions are small, and that the density fluctuations are related to the velocity 
fluctuations according to 


<p’> 

P 


(Y-l) M a 2 


<u' > 

U 


These assumptions appear to hold even in severely perturbed boundary layers 
(see, for example Dussauge and Gavigllo, 1981). In addition to the above, 
the value of the density-velocity correlation is required. The measurements 



-18 


i 




by Dussauge and Gaviglio suggest a constant value of 0.8 across the boundary 
layer, and this value was adopted for all data analysis. 

The inclined wire measures the mass-weighted shear stress (pu) 'v' . 

To determine the behavior of the turbulent shear stress T « pu'v', we again 
assumed that the "Strong Reynolds Analogy" could be applied. That is. 


t/t 


w 


p u'v* _ _2_ p U ^ 1 (pu) ' v* 

p w U T 2 Cf Pref^ref [l + ( Y-D M a^ P ^ 


where the subscript "w" indicates at the wall. In the case of the undisturbed, 
upstream boundary layer (see Fig. 17), the results agree closely with Sandborn's 
"best fit", which does not prove that the data are accurate, but nevertheless 
gives further confidence in our method. 

Before presenting the turbulence measurements, it is useful to consider 
some aspects of the mean flow behavior. 

The static pressure and skin frictiovj coefficient distributions are 
shown in Figs. 18 and 19, respectively. From these figures, it is clear 
that the flow tends towards separation as the corner angle increases. In 
the 20° case, the flow has actually separated, with a separated zone approxi- 
mately one-half 5 long. 

The mean velocity profiles generally demonstrate a quick recovery 
downstream of the corner. For all three cases, Settles, et al. (1979) ob- 
served that the profiles rapidly fill out and approach their equilibrium 
shape by the furthest downstream station. 

3.2.1 Normal wire results 

When we consider the turbulence behavior, however, it quickly becomes 
clear that the boundary layer is far from equilibrium* even at the furthest 
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point downstream. Figures 20 and 21 show the mass-flow fluctuation intensity 
<(pu)’> and velocity fluctuation intensity <u’>. Bote that the intensities 
are non-dimensionalized by the incoming freestream values (pU) re f and U re f • 

The figures therefore show the behavior of the absolute fluctuation levels. 

The most obvious feature of the turbulence behavior is the dramatic 
amplification that occurs as the boundary layer passes through the interaction 
region. To understand this amplification more fully, it is useful to consider 
the evolution of the turbulence intensities along selected streamlines. In 
this way, the amplification may be related to the terms which appear in the 
Reynolds stress transport equations (see, for example, Hayakawa, et al., 1983). 

Three streamlines, corresponding to upstream boundary layer locations 
of y/<$ = 0.2, 0.4 and 0.6 were selected for this purpose (see Fig. °.2) . For 

each of these streamlines, the evolution of u'2, normalized by its upstream 
value is given in Fig. 23. 

Figure 23 shows that these flows may be conveniently divided into two 
regions: an "interaction zone," which loosely corresponds to the region where 
severe pressure gradients exist, and a "recovery zone, " which is the region 
downstream of the interaction zone. Within the interaction zone, u'2 dis- 
plays a rapid increase along each streamline. The magnitude of the increase, 
however, is a function of both the initial position within the boundary layer 
and the shock strength; the largest increase occurs for the largest corner 
angle. In contrast, the behavior within the recovery zone does not display 
such a simple trend. Near the wall, the turbulence Intensity decays quite ra- 
pidly (Fig. 23a), whereas further away the intensities are still increasing, even 
at the furthest downstream station (Fig. 23b). This Behavior clearly reflects 
the time dependent response of the turbulent motions,. As expected, the. larger 
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eddies in the outer part of the boundary layer obviously respond more slowly 
than the motions near the wall. 

3.2.2 Inclined wire results 

The mass weighted Reynolds shear stress (pu) 'v'/(pU^) re f is shown 
in Fig. 24. Before discussing these results, the accuracy and limitation 
of the hot-wire technique should be considered. 

It is becoming increasingly clear that one of the major limitations 
on using inclined hot-wires in low Mach number supersonic flows is the 
requirement that the instantaneous normal Mach number should be greater 
than 1.2 (Smits and Muck, 1984) . Within the transonic range, the hot-wire 
behavior changes drastically, and a sharp fall-off occurs in the inferred 
turbulence intensity. For this reason, results for which the normal Mach 
number is below 1.2 are expected to be in error and should be used cautiously 
although it may be possible to use the shear stress at the wall as a guide 
for interpreting the results near the surface. 

Consider now the results shown in Fig. 24. ’Re behavior of (pu)'v' is 
obviously rather different from that displayed by For example, in 

passing through the interaction region, (pu)'v' for the 16° corner increases 
by as much as 16 times, whereas ( u)'2 increases by less than 10 times over 
the same distance. For the 20° corner, the maximum level at the last station 
is more than 20 times its upstream value. This may le comparer, to (pu)'^ 
which increased by about 16 times over the same distance. 

The behavior of either the Reynolds stress pA', or the kinematic 
shear stress u'v' may be deduced from the results foar (pu)'v' by using the 
Strong Reynolds Analogy. Here we have chosen to present only the results 


for the kinematic stress u'v'. 
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Figure 25 shows u'v'/U ^ plotted against y/& Q for the different 

streamwise stations. The behavior of u'v'/U re £ is qualitatively similar to 

? — 2 

that of u ' z / u re £* Quantitatively, however, the amplification of u'v' is not 

2 

as great as that observed for u’ ; this suggests that the turbulence structure 
is significantly changed through the interaction. 

For instance, the ratio - u'v'/u'^ can be taken as a structure para- 
meter. In the undisturbed boundary layer at x = -51 mm, this ratio has a 
value of about 0.25 (at y/v = 0.6), which agrees well with the value commonly 
quoted for incompressible boundary layers (see Townsend, 1976, p.107). Through 
the interaction zone of the 16 corner, the ratio u’v'/u'^ increases to about 
0.32, although it must be stated that the measurements of u'v' in this region 
are probably not too reliable. Downstream in the recovery region, however, 
this ratio decreases significantly, and it reaches a value of 0.16 't 0.18 at 
the furthest downstream station (x = 140 mm). A similar trend was observed 
by Smits et al. (1979) in a subsonic boundary layer subjected to an impulse 
in concave curvature. It is interesting to note that in the present experi- 
ment, where both curvature and compressibility effects are important, the 
turbulence appears to respond in a fashion similar to that observed in. an 
incompressible curved flow. 

In passing the interaction region of the 20° corner, however, 
u'v'/u '2 drops by about 25%, and remains low, although it appears to recover 
slightly at the last two stations. This result contrasts distinctly with 
the behavior observed in the 8° and 16° corner studies: there the ratio 

Increased sharply through the interaction before relaxing somewhat downstream. 
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3 . 2.3 Final discussion of compression corner flows 

All turbulence quantities (pu) , u' , (pu) 'v' and u'v' show a 

qualitatively similar trend; all increase steeply on encountering the shock 

and continue to increase over a distance which corresponds approximately to 

the region of non-zero pressure gradient. The relaxation process downstream 

appears to be rather slow, and at the furthest downstream station (about A 6 Q 

downstream of the corner), the turbulence is still far from equilibrium. 

Perhaps the most remarkable feature is that for the 20° comer, the dimension- 

— — ^ 

less turbulence structure parameter u'v'/u’' is significantly reduced in 
passing through the interaction, in contrast with behavior observed with 8° 
and 16° corner flows. 

To interpret these results, we can begin by considering the theoretical 
work of Zang et al. (1982) and Anyiwo and Bushnell (1982). Using the two- 
dimensional Euler equations, these authors showed that several possible turbu- 
lence amplification or generation mechanisms may occur during an interaction 
between a plane shock wave and an incident turbulence field. These mechanisms 
include (1) direct amplification, (2) "generation" of turbulence from incident 
acoustic and entropy fluctuations, and (3) "focussing" caused by distortions 
of the shock front. 

When a shock wave interacts with a turbulent boundary layer, however, 
the flow field is rather different from the ideal case considered by these 
authors. The shock wave within the boundary layer curves and is followed 
by a system of instantaneous compression waves (see, for instance. Fig. 3). 

All these waves are unsteady and appear to be in constant, low frequency 
motion (Dolling and Murphy, 1982). As Zang et al . (1982) point out, the 
direct conversion of mean flow energy into turbulence by shock oscillation 
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can serve as a powerful turbulence amplification mechanism. 

For both the 8° and 16° corner flows, the effects of dilatation 
and curvature appeared to be sufficient to explain the measured amplification. 
Direct shock effects and shock oscillation did not seem to be overly important. 
This conclusion was supported principally by the behavior of the structure 
parameter u'v'/u . For example, Debieve et al. used an analysis based on 
rapid distortion theory to show that this parameter should increase through 
a shock, when the shock behaves as a stationary discontinuity. In addition, 
Smits, et al. (1979) found that u'v'/u'^ increased when a subsonic boundary 
layer was subjected to a short region of concave curvature. In the 20° 
corner flow, however, this parameter decreased through the interaction. To 
explain this result, we tentatively suggest that shock oscillation may be im- 
portant in this flow, in addition to the effects of dilatation and curvature. 
Shock oscillation is likely to become more important at these high corner 
angles because the strength of the shock increases and eventually produces 
separation. If the shock movement is essentially random, we expect that the 
mean flow energy is transferred more to the normal stresses rather than to 
the shear stresses. This explanation suggests that the motion of the shock 
wave generates significant "inactive" motions such as those discussed by 
Bradshaw (1967) in relation to highly retarded subsonic boundary layers. 

3.3 Isentropic Compression Studies 

We examined the evolution of the mean flow and the turbulent stresses 
in a supersonic boundary layer experiencing the effect of bulk compression 
and streamline curvature. This study was prompted in part by the investiga- 
tions of shock-wave/boundary layer interactions generated by two-dimensional 
compression comers (see Section 3.2). In that work, it was observed that 


s 


the boundary layer parameters suffered considerable distortion; in particular, 
the turbulent stresses were dramatically amplified. Four separate mechanisms 
for turbulence amplification were identified: "direct" amplification by the 

shock wave oscillation; the destabilizing effect of compression downstream 
of the shock, and the destabilizing effect of concave curvature. 

The present investigation was designed to isolate the effects of 
compression and curvature. Instead of turning the flow suddenly, as in a 
compression corner, the turning was accomplished more gradually, through a 
short region of concaveiy curved wall. By spreading the pressure rise over 
several boundary layer thicknesses, the shock wave forms outside the boundary 
layer and it has no direct effect on the boundary layer behavior. 

The incoming boundary layer and upstream flow conditions were 
identical to those used in the compression ramp studies. The upstream 
boundary layer was allowed to develop fully under a nominally zero pressure 
gradient on the wind tunnel floor before entering a short region of surface 
curvature. Just upstream of the curved wall, the boundary layer thickness 
<5 0 was approximately 26 mm. The total turning angle was fixed at 8°. Two 
different constant radii curvatures were investigated: the first had a 

radius of 254 mm, the second a radius of 1270 mm. These correspond to ratios 
of $ 0 /R of approximately 0.1 and .02 respectively. Both curvatures were 
followed by a short (153 mm) recovery region, allowing us to study the 
initial relaxation behavior of the boundary layer (see Fig. 26). 

The mean flow behavior for both models was previously reported by 
Taylor and Smits (1984) and Taylor (1984). Briefly, the velocity profiles 
displayed a "dip" below the log-law, suggesting an increase in the turbulence 
length scale, and no evidence was found for the presence of longitudinal 



roll-cells, as might be expected to occur in concavely curved shear layers. 

In both these respects, the flow was similar to the incompressible flow 
studies by Smits et al. (1979), who investigated the response of a -boundary 
layer as it experienced the combined effects of concave curvature and lateral 
divergence. In the present flow, the boundary layer experiences the combined 
effects of concave curvature and bulk compression, and the analogy between 
divergence and compression, first suggested by Green in Bradshaw (1973), was 
strengthened considerably by these mean flow observations. 

Measurements of the turbulence behavior were also taken. Normal and 
inclined hot-wires were used, and the measurements describe the behavior of 
the longitudinal mass-flow fluctuations (pu) ' and the mass weighted shear 
stress (pu)'v'. The kinematic stresses u'^ and u'v' were deduced using 
Morkovin’s "Strong Reynolds Analogy." Wherever possible, the results for 
the 8° compression comer were compared with those for Models I and II. The 
incoming boundary layer and the upstream freestream conditions were identical 
in all three experiments. The static pressure and skin friction distributions 
are compared in Figs. 27 and 28, respectively, and the figures demonstrate 
that the experiments cover a wide range of stress gradients, that is, a wide 
range of strain rates, although the overall pressure rise and turning' angle 
are the sane in each case. Hence, these experiments investigate, for a given 
perturbation strength, the effect of varying perturbation rate. 

The variation of the longitudinal component of the turbulence intensity 
is shown in Fig. 29. The upstream reference conditions were used as non- 
dimensionalizing variables throughout, and hence the results show the behavior 
of the absolute turbulence levels. 
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The general behavior of u'^ is very similar in all cases. Initially, 
there is a rapid amplification of u'^ through the perturbation zone. This 
increase continues further downstream, and the peak amplification of u'^ is 
about three-fold in each case. In the outer part of the boundary layer, 
there is little sign of any recovery or relaxation but near the wall the 
intensities quickly fall, indicating that the relaxation process begins at 
the wall and propagates outward. 

The variation of the kinematic shear stress u'v'/U* . is shown in 

ref 

Fig. 30. Again, a considerable amplification occurs through the perturbation 
zone, continuing into the region downstream, and the relaxation appears to 
propagate outwards from the wall. However, the peak amplification levels 
for the three cases differ considerably; for the compression corner and 
Model I it is about four-fold, whereas for Model II it is only about two-fold. 

Hence, the present results show that structure parameters, such as 
u'v'/u'^, are clearly a function of the perturbation rate. 

As far as the turbulence behavior is concerned, the compression corner 
and Model I influence the turbulence in an almost identical manner. This 
result suggests that the perturbation in these two cases is sufficiently 
rapid to alter the turbulence in a manner which depends only on the overall 
changes that occur, not. on the path taken. The parameter describing this 
"change of state" is therefore the total strain, that is, the integral of 
the strain rate over the time it acts. For the curvature, this integral is 
equal to the total turning angle, and for the compression it is equal to 
(1/y) f.n(p 2 /pi) . It appears that the shock wave itself has no explicit effect 
on the turbulence but this conclusion can only hold if the shock is, relatively 
weak, such that the entropy loss is small, and separation with the concommitant 
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unsteadiness is avoided. This condition seems to be satisfied in the 8° 
compression corner experiment. 

The response of the turbulence to the perturbation produced by Model II 
is quite different to that seen in the other two experiments. For example, 
the stress ratio u'v'/u' z is affected by only a relatively small amount. 

Hence, despite the large amplification of the absolute turbulent stresses, 
the turbulence structure is not radically altered, implying that the perturba- 
tion is sufficiently slow for some redistribution processes to occur. In this 
case, the local strain rates are probably more useful than the total strain 
for describing the response of the turbulence. The results from this parti- 
cular experiment, therefore, appear to constituce an excellent test case for 
developing turbulence models. 

Close to the wall, we observed the beginning of a relaxation process 
in all three experiments. It is to be expected that the flow near the wall 
will attain equilibrium more quickly than the flow in the outer part of the 
layer; a measure of the large eddy time constant is the turbulent energy- 
divided by its rate of production, and this will vary approximately as 
l/(ol’/5y) (Bradshaw, 1973). This relaxing region grew in size as we pro- 
ceeded downstream, and its growth resembled that of a new boundary layer. 
Similar "internal layers" have been observed in boundary layers perturbed 
by sudden changes in surface roughness, surface curvature, and pressure 
gradient (Smits and Wood, 1985) and the similarity displayed by the propa- 
gation of the relaxation outward from the wall may be useful in a qualitative 
understanding for the flow behavior in the present experiments. 
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3.4 Summary and Conclusions of Experimental Work 

When taken as a whole, the experimental work leads to several interesting 
conclusions. For example, we have seen that turbulence levels are strongly 
amplified in a shock-wave boundary layer interaction. This amplification 
appears to be caused' by direct, virtually inviscid amplification across the 
shock, followed by the combined effects of adverse pressure gradient, compres- 
sive extra strain-rates and concave curvature. When the shock strength is 
relatively low (that is, no separation occurs) it seems that the important 
parameter is the overall pressure rise rather than the presence of a shock- 
wave. This was demonstrated by the results from the isentropic compression 
studies. When the shock strength increases, however, shock-wave oscillation 
becomes an important amplification mechanism. Here, mean flow energy is 
directly transferred into unsteady turbulent motions. This process is 
apparently random, and therefore contributes more to the random motions 
which increase the total turbulent energy than to the organized motions- 
associated with the shear stresses. This hypothesis explains why the struc- 
ture parameter u’v'/u'^ increases' through the interaction for the 8° and 16° 
compression corners, and decreases for the 20° compression corner; when. the 
flow separates at higher angles, shock-wave oscillation becomes important. 
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4. ANALYTICAL AND COMPUTATIONAL WORK 


Throughout the period of our previous funding, we have maintained 
close contact with several computational groups, especially those groups 
headed by Prof. D. D. Knight at Rutgers University, and Dr. C. C. Horstman 
at NASA-Ames. In addition, our contacts include R. E. Melnik (Grumman Aero- 
space Corporation), Ha Minh Hieu (l'Institut National Polytechnique de Toulouse), 
D. Degani (Technion), J.-P. Bonnet (Poitiers) and H. H. Legner (Physical Sciences 
Inc.). All of these groups have expressed interest in trying to calculate the 
compression corner flows, although they generally use rather conventional tur- 
bulence models. 

What is perhaps more promising is the approach suggested by Dussauge 
and Gaviglio (1981), Debieve et.al. (1982) and Anyiwo and Bushnell (1982). 

In this work, the concept of "sudden distortion" is applied to compressible 
flows experiencing a short region of intense pressure gradient, including 
the interaction with a shock-wave. When the perturbation is sudden, that is, 
the perturbation occurs over a time which is considerably less than the 
response time of the turbulent motions, then the distortion may be essentially 
inviscid. The turbulent motions are then affected only through the distortion 
of the mean field, and nonlinear effects caused by turbulence-turbulence 
interaction can be neglected. This approach appears to have considerable 
promise in calculating the turbulence behavior observed in our experimental 
work. We have been fortunate enough to have had Dr. J.-P. Dussauge from 
I.M.S.T. in Marseille as’a Visiting Research Fellow this year, and we are 
working closely with him to complete the calculation of our compression corner 
flows as well as the "isentropic" compressions. Preliminary work suggests that 
this essentially inviscid analysis might supplement more conventional schemes, 
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5. MISCELLANEOUS ADDITIONAL WORK 


Earlier this year. Professor H. Femholz from the University of Berlin 

! 

j spent six weeks in the Laboratory as a Visiting Research Fellow. During this 

I 

i period, we had the opportunity to interact closely with him, and we intend to 
maintain these discussions in the future. In particular, we are .raking our 
data available to him for future publication in AGARDograph form. 

In another development this year. Professor A. Smits of the Gas 
Dynamics Laboratory was invited to write a review article for Annual Reviews 
of Fluid Mechanics. Entitled "The response of turbulent boundary layers to 
a sudden perturbation," the article is co-authored with Dr. D. H. Wood from 
Newcastle University (Newcastle, Australia), and is due to appear early in 
1985. 
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6. PROPOSALS FOR FUTURE WORK 

The conclusions and observations from our previous work lead us to 
propose the following. Firstly, we are interested in testing our current 
understanding of shock-wave boundary layer interactions. Although we have 
identified some of the important physical mechanisms, we are still far short 
of a quantitative understanding. Further theoretical and experimental work 
is required to produce a useful turbulence model. In the experiments per- 
formed thus far, many different effects were acting simultaneously. To 
understand these effects fully, we need to study them in isolation, and our 
first priority is to investigate relatively simple geometries such as a flat 
plate boundary layer subjected to a severe adverse pressure gradient. In 
this geometry, the effects of the shock, as well as effects of the curvature 
would be eliminated. Preliminary work has already commenced, and we have 
designed models to reproduce on a flat plate the pressure gradients observed 
in the isentropic compression corners. By comparison between these sets of 
experiments, we will be able to distinguish between the effects of compression 
and curvature. In a complementary study, we propose to complete the isentropic 
compression work by investigating a flow with a total turning angle of 10° and 
6 0 /R = 0.02. This work will provide important information on the range of 
applicability of Rapid Distortion Theory. The model has already been constructed, 
and the mean flow studies have been completed (publication [22], Appendix A), 
although the turbulence measurements have not yet commenced. 

In addition, we are interested in studying the time-dependent behavior 
of turbulent motions in a supersonic boundary layer. We feel that this work 
is essential to a proper understanding of turbulent transport properties, and 
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it will represent a significant departure from our previous work which used 
the more conventional approach based on Reynolds-averaged quantities. 

We propose to begin with a detailed study of the 24° corner at Mach 
2.9. Here, a substantial separated region exists, and shock-wave oscillation 
is expected to be particularly important. Currently, Drs. K. C. Muck and 
J.-P. Dussaugt are measuring the correlation between wall pressure and velo- 
city fluctuations. These measurements will provide information regarding 
the "inactive" motions, and the data should test our suggestion that shock- 
wave oscillation amplifies the random motions more than the organized motions. 
In this study, we have also used an array of four wall pressure transducers 
to study space-time correlations of the pressure signal and the distribution 
of phase velocities. Preliminary results indicate that the shock-wave oscil- 
lation displays both a bulk motion backwards and forwards as well as a spanwise 
wrinkling effect. We propose to complete the analysis of existing data 
and to extend the experimental work. In this work, we propose to begin 
with a set of conditionally sampled measurements wehre the shock wave posi- 
tion is used as a trigger for the conditional sampling. In this way, the 
instantaneous conditions before and after the shock wave can be determined. 

We also propose active control of the shock wave position by applying a 
perturbation to the separation bubble. By controlling the shock wave motion, 
the effect of shock wave oscillation may be more clearly defined. To study 
the instantaneous wall conditions, an array of thin hot -film gauges will be 
mounted flush with the wall. These gauges will provide instantaneous wall 
shear information, and will enable us to discern the connection between the 
shock wave motion and the behavior of the separation bubble. 

On a more fundamental level, we propose a detailed study of the large- 
scale motions in a flat plate sero pressure-gradient boundary layer. A great 
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deal of work has been performed to identify the characteristic motions, 
sometimes called coherent structures, which occur in turbulent shear flows, 
but our current understanding of these large-scale motions is mainly derived 
from subsonic flow studies at relatively low Reynolds number (sec, for example, 
Brown G Thomas, 1977). We propose to use multiple hot-wires, combined with 
arrays of flush-mounted hot films, to similarly characterize the nature of a 
compressible , high Reynolds number boundary layer. We believe that this 
study would be the first of its kind. Some work has already commenced in this 
area, although the measurements thus far have been confined to mapping out the 
overall nature of the boundary layer. By measuring uv and uw in the spanwise 
and streamwise directions, we have established that the boundary layer is 
acceptably uniform and two-dimensional. Preliminary measurements of the wall 
pressure-velocity correlations are currently being analyzed. 

In another major effort, we propose to examine the relaxation behavior 
of the boundary layer far downstream of a severe perturbation such as an 
interaction with a shock wave. This behavior is of great practical interest. 
For example, as the flow enters an inlet, the boundary layers interact with 
the entry shock system and then relax before entering the compressor stage. 
Thus, it is the relaxation behavior which really governs the engine inlet 
condition. The relaxation process is highly nonlinear, however, and very 
difficult to calculate. It is also experimentally difficult, and no measure- 
ments are currently available far downstream. In our present facility, such 
measurements are not possible but our new facility (funded by a recent DOD 
equipment grant, see Appendix C), will allow this work to proceed. This 
facility features a flexible diffuser, which by careful adjustment will 
create a relatively long disturbance-free region downstream of the interaction 
(> 50 6 0 ) . 
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We propose to begin with a study of the relaxation behavior downstream 
of an 8° compression corner. In terms of turbulence modelling, our ultimate 
aim is to combine a rapid distortion analysis for the interaction region with 
a suitable transport model (such as a k-c model, for example) to describe the 
relaxation process. 

As support for the work mentioned above, we intend to 
modify existing instrumentation and develop and test a variety of new 
instruments. The proposed work in hot-wire anemometry, and in the measure- 
ment of wall pressure fluctuations has already been mentioned. In particular, 
we wish to develop a hot-wire rake capable of taking velocity data simultaneously 
at a number of points within the boundary layer. We have also alluded to flush- 
mounted thin-film gages for the measurement of instantaneous wall shear stress. 
These gages were suggested by Prof. Nosenchuck of our Department, and have the 

potential to measure the mean and fluctuating skin-friction, + t w \ 

Several other techniques are being developed within the Laboratory, 
with particular emphasis on new and improved flow visualization methods. 

Our long-range goal is to visualize the whole flowfield and its behavior as 
a function of time, with sufficient resolution to examine detailed areas of 
the flow. We have identified many possible techniques, including sharp- 
focus schlieren, multiple, schlieren, tracing "hot-spo^s" generated by spark 
discharges or high-energy laser pulses, multiple high-speed imaging techniques 
and local vapor screen methods. 

We are also continuing our development of the Resonant Doppler 
Velocimeter (RDV) as a flow visualization tool. Here, the fluorescence of 
an atomic or molecular seeding species is used to highlight areas of high 
or low velocity, temperature and pressure (see Zimmerman and Miles 1983 and 
Zimmerman et al. 1983). 



In the long term, the RDV techniques will allow quantitative analysis 
of mean and fluctuating pressure, velocity and temperature within the flow 
field (Cheng et al. 1983). In the short term, however, we feel that another 
technique is required to supplement (and cross-check) the hot-wire measurements. 
We expect, therefore, to begin laser-Doppler measurements within the next three 
years. This will require considerable development, plus capital investment, 
and we hope to find support for this work. 

One of our primary purposes in pursuing this experimental work is to 
develop a quantitative model for the turbulence behavior in compressible 
boundary layers. Ke propose to begin this task by (a) extending the work 
performed by Bradshaw on small extra strain rates and impulsively applied 
extra strain rates, and (b) modifying the rapid distortion approach suggested 
by workers at the Institut de Mechanicque de la Turbulence at Marseille. 

In addition, the effect of shock wave oscillation will be examined by extending 
the work of Anyiwo and Bushnell (1982) and the work of Debieve et.al. (1982). 

We feel that these inviscid approaches to suddenly perturbed flows hold great 
promise, and we anticipate making substantial progress in this area. 
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Figure 1. Comparison between experiment (Settles et. al 1978) and compu- 
tation (Visbal and Knight 1983) for compression corner flows. 
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Sketch of current probe design. Wire is tungsten, 2 to 5 urn 
diameter, 0.8 cm long. 
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3 Detail of probe design, allowing approximate ehock locations, 
(Major chocks only.) 
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Fig. 5 a-c. « rms mass-flow fluctuation inlenstty for zero pressure gradient 2.9 turbulent boundary layer (R„ - '"’.MIO. d*» 2f mm. 
C, “0.0011). Results are for three different w ires with overheat ratio ^ 1.0 b nns mass-flow- fluctuation intensity for zero p:e>-u:c 
gradient i/„=2.9 turbulent boundary layer (/?, = 77,bOO. <5= 26 mm. C, = O.lfflit 1 ). Results are for a single wire, with varying overheat 
ratio, c A comparison between the inferred rms velocity fluctuation intensity (deduced from the results given in Fig. I0.il and some 
representative data reproduced from Dussauge and Gaviglio ( I9S1 ) 
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A sketch of a typical inclined wire probe. The tungsten 
wire was 5pm diameter. 
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p]g |Q Measurements of (pu)'v' and pu'v' in a zero pressure gradient 

turbulent boundary layer at Mach 2.9: a ; 30° wire, A ,45° wire 
Shaded area corresponds to Sandborn's ''best estimate" for the 
supersonic shear stress distribution (Sandborn, 1 9 74 ) . Error 
bars are only shown for the 30° wire, (t/t = (p u* v' )/p u 2) 
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Fig. 11 Calculated results for the directional sensitivity of 
nominally 30° and 45° wires (£ ■ 1 mm, d = 5pm) : 

, infinitely long wire (equation 2); 

, straight wire with end conduction; 

l t bowed wire with end conduction (6 = ±10) 

(a) subsonic flow (Re = 4, p = 1.0), (b) supersonic 
flow (Re *= 100, p = 2.0). 
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Fig. 16 Sketch of the 20 x 20 cm High Reynolds 
Number Channel 
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Fig. 17. Profiles of the Reynolds shear 
stress normalized by the surface shear 
stress, t/t„. for different wires at the 
upstream station (x = - 50.8 mm). 
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